Separability criteria for continuous variable systems 
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A general separability condition on the second moment (covariance matrix) for continuous variable 
two-party systems is derived by an analysis analogous to the derivation of the Kennard's uncertainty 
relation without referring to the non-negativity of the partially transposed density matrix. This 
separability criterion is generally more stringent than that used by Simon which is based on the 
non-negativity of partially transposed density matrix, and thus this criterion may be useful in the 
analysis of general continuous two-party systems. Another separability criterion used by Duan et al. 
is shown to be generally weaker than that of Simon. We thus have a hierarchy of separability criteria, 
but all these criteria when combined with suitable squeezing become equivalent at the boundary of 
the P-representation condition and thus turned out to be sufficient to analyze the separability of 
two-party Gaussian systems. 

PACS numbers: 



I. INTRODUCTION 

The entanglement [l[ is an intriguing property of quan- 
tum mechanics, but a general quantitative criterion of en- 
tanglement appears to be missing at this moment. The 
negativity or non-negativity of the partially transposed 
density matrix was proposed as a quantitative means to 
analyze the entanglement by Peres [3] • But this condition 
is known to work only for simple systems such as a two- 
spin system [3| . In view of this situation, it is remarkable 
that the necessary and sufficient conditions exist for the 
separability of two-party Gaussian systems fj, Q. To 
be precise, the "two-party systems" in the present pa- 
per are used for the systems with one freedom (or one- 
mode) in each party, since it is known that the analysis of 
the separability criterion for the two-party systems with 
more than one freedom in each party is more involved 
@. The criterion used by Simon Q is based on a gen- 
eralization of the non-negativity condition of the partial 
transposed density matrix in the manner of Peres. On the 
other hand, the criterion used by Duan et al. [5| is based 
on a set of relations for EPR-type operators. The non- 
negativity condition was further analyzed by Shchukin 
and Vogel @, H[. This separability and related issues 
have been discussed in the past by various authors, for 
example, in @, EE El El El EI El EE E3 ■ The present 
status of the quantum separability problem of two-party 

See also a review of 
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Gaussian states is reviewed in 
continuous variable systems in 

We here start with an elementary analysis of the 
Heisenber g u ncertainty relation in the manner of Ken- 
nard Mi \M and we derive a separability condition, 
which is generally more stringent than the separability 
condition given by Simon Q , without referring to the neg- 
ativity of the partially transposed density matrix. This 
implies that our criterion is more efficient as a condition 
for separability than the condition used by Simon, which 
is based on the non-negativity of the partially trans- 
posed density matrix, although neither criterion provides 
the sufficient separability condition for two-party systems 



in general. Our analysis is complementary to that of 
Shchukin and Vogel 0, Q ; they analyze higher moments 
also in the framework of the non-negativity of the par- 
tially transposed density matrix, while we specialize in 
the second moment and search for a more stringent con- 
dition. We also show that another class of separability 
conditions used by Duan et al. [|[ are generally weaker 
than the condition of Simon. Despite these differences, 
all these criterions when combined with suitable squeez- 
ing give rise to the necessary and sufficient condition for 
the separability of two-party Gaussian systems. This fact 
has been shown elsewhere by using the analytic solutions 
of squeezing parameters [22j which establish the equiva- 
lence of the P-representation condition for Gaussian sys- 
tems with the Sp(2, R) ® Sp(2, R) invariant separability 
condition. Some of the technical aspects involved in the 
analysis in [22| are also clarified in the present paper. 



II. ENTANGLEMENT AND KENNARD'S 
RELATION 

A. Kennard's relation 

We consider a two-party system with one freedom in 
each party (or a two-particle system in one-dimensional 
space) described by canonical variables (<7i,pi) and 
(<}2,P2)- We define 

X(d, f) = dtqi + d 2 pi + fiq 2 + f 2 p 2 , 
X(g,h) = gxqx + g 2 pi + + h 2 p2 (1) 

where all the coefficients 

d T = (d 1 ,d 2 ), / T = (/ 1; / 2 ), g T = ( gi ,g 2 ), 

h T = {h 1 ,h 2 ) (2) 

are real numbers. The use of the operators in (1) is known 
to be convenient to define the general uncertainty rela- 
tions. In fact, the conventional form of the Kennard's 
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relation is derived as a special case, as is shown in (14) 
below. We deal with a mixed state p — ^ fc -Pfc|V'fc)(V'fc| 
with P k > and E fe -Pfe = !• 

The Kennard's relation for a specific pure state | } is 
written for any choice of parameters d to h as 

((X(d, /) - ,/)) fe ) 2 ) fc + ((j£( 5 , ft) - (X(g, h)) k f) k 
>\d T Jg + .f T Jh\ (3) 

where the 2x2 matrix is given by 

1 



J = 



-1 



(4) 



and (X(d,f)) k = {ipk\X(d, f)\ipk), for example. The 
relation (3) is derived from (ipk\vtf l^k) > and 
(ipkWvliPk) > for 

f, = X(d, /) - f)) k + i(X(g, h) - (X(g, h)) k ), (5) 

and the right-hand side of (3) stands for the commutator 

\[X(d,f) - (X(dJ)) k J(g,h) - (X(g,h)) k }\ 
= \d T Jg + f T Jh\. 

We thus conclude by taking the weighted sum of (3) 
J2Pk((X(d,f)-(X(d,f))k) 2 )k 

k 

+ J2Pk((X(g,h)-(X(g,h))k) 2 )k 

k 

>\d T Jg + f T Jh\ (6) 

On the other hand, it is customary to deal with 

((AX(d,f)f)=Tr{(AX(d,f)fp} (7) 

where AX(d, /) = X{d, f) - (X(d, /)) with 

(X(d,f))=Tr{X(d,f)p} (8) 

in the analysis of uncertainty relations for a mixed state. 
We note the relation 

((AX(dJ)) 2 ) = ((x(d,f)-(X(d,f))) 2 ) 

= J2 p k((x(d,f)-(X(d,f))) 2 ) k 
k 

= J2 P k({X(d,f)-(X(dJ))k 

k 

+ (X(d,f))k-(X(dJ))} 2 )k 

J2Pk[((x(dj)-(x(d,f)) k ) 2 )k 

([X{dJ))k-{X{dJ)))\ (9) 



From (6) and (9) we thus conclude the basic relation 

((AX(d,f)) 2 ) + ((AX(g,h)) 2 ) 

= Y,Pk[{{x(d,f)-{X(d,f))k) 2 )k 

k 

((X(d,f))k-(X(d,f))) 2 } 

-J2Pk[((x(g,h)-(X(g,h))k) 2 )k 

k 

((X{g,h))k-(X(g,h))) 2 ] 

>J2Pk((X(d,f))k-(X(d,f))) 2 

k 

-J2Pk((X(g,h)) k ~(X(g,h))y 



+ 



+ 



+ 



+ \d T Jg + .f T Jh\. 



(10) 



This relation is more precise than the customary form of 
the uncertainty relation for a mixed state 

((AX(d, J)) 2 ) + ((AX(g, h)) 2 ) > \d T Jg + f Jh\ (11) 

which is obtained from Tr{fj'(fj')^ p} > and 
Tr{(f]'yf]'p} > with 

ff = AX(d,f) + iAX(g,h), 

and the right-hand side of (11) standing for the commu- 
tator |[AX(ci, /), AX(g, h)]\. Note that one may choose 

(X(d,f)) - Ek P k(X(d,f))k - 0, but (X(d,f)) k for 
each component state with 

(X(d,f)) k = J d qi dq 2 Vk(qi,q2)X(d,f)M<li,Q2) (12) 

does not vanish in general in (10). 

In passing, the Kennard's relation for a general pure 
state (3) is also written as 

x 2 ((x{dJ)-(X(dJ)) k ) 2 )k (13) 
+((x(g, h) - (X(g, h))k) 2 )k > x\d T Jg + f T Jh\ 

by replacing d — > xd, f — > xf for any real positive x 
(and this relation holds for negative x also) , and thus the 
discriminant gives the conventional form of Kennard's 
relation 



fe 
+ 



((x{d,f)-(Jt(d,f)) k ) ) k ^{(x{g,h)-{X(9,h))k) )k 
>\\d T Jg + f T Jh\ (14) 

if the parameters d to h are suitably chosen. 
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B. Separability condition 

For a separable pure state ipk(qi,Q2) = </ , fc(<Zi)v 3 fc( ( 72), 
we have 



where = £ M - (^) with = J2k p k(^Pk\^\ipk) ■ 
We here recall the definition of the density matrix p = 
Z)fe Pk\ipk)(ipk\ with P fe > and £) fc Pk = 1- The corre- 
lation matrix can be written in the form 



((i(d,/)-(i(d,/)) t ) 2 ) fc 

= ((rfigi + d 2 pi - (digi + d 2 pi} k ) 2 ) k 

+ <(/l«2 + /2P2 - (/l«2 + f2P2)kf)k (15) 

since the cross terms vanish, and similarly for 
{(^X(g,h) — (X(g,h))kj )k- Wc thus have for a sepa- 
rable pure state 

((X(dj) - (X(dJ)) k ) 2 ) k + ((x(<7,/i) - (X{g,h))k) 2 )k 

= {(x{d,o)-{x{d,o)) k ) 2 ) k 

+ ((x(g,0)-(X(g,0))k) 2 )k 

+((x(o,f)-(x(o,f)) k ) 2 ) k 

+ (( y X(0,h)-(X(0,h)) k y) k 
> \d T Jg\ + \f T Jh\ 



V 



A C 
C T B 



(19) 



where A and B are 2x2 real symmetric matrices and C 
is a 2 x 2 real matrix. We also define 

V = {%„), = J2 Pk(&L)k(AL)k (20) 



and 



V = 



A C 
C T B 



(21) 



(16) 



which holds for any choice of parameters d to h. Here we 
used (3) for / = h = or d = g = 0. The equality sign 
holds only for 

[(diqi + d 2 p\ - (diqi + d 2 pi) k ) 

+i (gm + gin - (gm + 92Pi)k)]<t>k(qi) = o, 
Khfc + hvi - (fm + hp2tk) 

+i {hiq 2 + h 2 p 2 - (hiq 2 + h 2 p 2 ) k )]ip k (q 2 ) = 0, 

for a suitable choice of parameters d to h. 

We finally conclude from the first expression in (10) 
when combined with (16) for any separable density ma- 
trix 

((Ai(rf,/)) 2 ) + ((AX(9,/ l )) 2 ) 
Y,Pk{((X(d,f)) k -(X(d,f))y 



> 

k 

+ ({X(g,h)) k -{X(g,h))y 
+ \d T Jg\ + \f T Jh\ 



where A and B are 2x2 real symmetric matrices and C 
is a 2 x 2 real matrix. Both of V and V are non-negative. 

The basic relation (17) for separable states, which is a 
necessary condition for separability, is then written as 

d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 

> d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 
+ \d T Jg\ + \f T Jh\. (22) 

while the Kcnnard's relation for general states (10) is 
written as 

d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 

> d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 
+ \d T Jg + f T Jh\. (23) 

Note the difference between |(i T Jg\ + \f T Jh\ and \d T Jg + 
f T Jh\. The appearance of the non- negative matrix V 
on the right-hand sides of these relations is the novel as- 
pect of our formulation. The anti-symmetric commutator 
parts in 

(A£>i) = i({A4,Ai}) + i([A4,Ai]), (24) 

which may be added to A and B in (19), do not con- 
tribute to (22) since the parameters d to h are all real. 

Under the Si <g> S 2 e Sp(2, R) <g> Sp(2, R) transforma- 
tions of (<7i,pi) and (q 2 ,p 2 ), respectively, we have 



(17) 



which holds for any choice of parameters d to h. This 
relation (17) is our basic necessary condition for the sep- 
arability of general two-party systems. 

We next define the variables = {qi,pi, q 2 ,p 2 ) and 
the 4x4 correlation matrix V by 



A 
A 



B 
B 



S 2 BS 2 , C 
S 2 BS 2 , C 



S±CS 2 
SxCST (25) 



V = (V), 



(18) 



V = i(A£>i + AiA^) = 1({A4, Ai}) 



which is equivalent to the transformation 

d^Sfd, f^Slf, g^Slg, h^S^h (26) 

in (22) if one recalls J = SxJSf and J = S 2 JS% 
; the inequality (22), which is valid for any choice of 
parameters d to h, holds after the transformation (26) 
and in this sense (22) is invariant under the above 
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Sp(2,R) <g> Sp(2,R). To be precise, we do not use any 
property of the wave function under Sp(2, R) <g> Sp(2, R), 
and thus our Sp(2, R)®Sp(2, R) transformation is rather 
defined by (25) for given constant matrices A, B and C . 

The difference between the two expressions in (22) and 
(23) appears when one replaces B and C by Sj BS3 and 
CS3 (and also B and C by S3BS3 and CS3), respec- 
tively, with 



1 

-f 



(27) 



One can undo the replacements in the expression (22) by 
transformations / — ► S3/ and /i — > S3/1, whereas it leads 
to |d T J<7 — f T Jh\ in the relation (23). The separabil- 
ity condition thus demands that the Kennard's relation, 
when written in the form (22), should hold both for the 
original system and for the system with the replacements 
of B and C by S3BS3 and CS3, respectively, which may 
a priori correspond to an unphysical density matrix for 
inseparable systems. By using S3, one can adjust the 
signature of detC at one's will [J] in the separability con- 
dition (22). 

It is also useful to consider the separability conditions 
weaker than the separability condition (22) by imposing 
subsidiary conditions g — J T d and h = J T f or g = J T d 
and h = —J T f, respectively, 

d T Ad + f T Bf + 2d T Cf + d T JAJ T d 
+f T JBJ T f ± 2d T JCJ T f 
> d T Ad + f T Bf + 2d T Cf + d T JAJ T d 
+f T JBJ T f ± 2d T JCJ T f 
+ (d T d + f T f) 

which is also written as 

A + JAJ T C ± JCJ T 
C T ±JC T J T B + JBJ T 

A + JAJ T C ± JCJ T 
C T ±JC T J T B + JBJ T 



(28) 



> 



(29) 



For general inseparable states in (23), we have only the 
first condition in (28) if one wants to keep (d T d + f T f) 
on the right-hand side in this form. The basic Sp(2, R) ® 
Sp(2, R) invariance of uncertainty relations (22) and (23) 
is lost in these weaker conditions (28) and (29) with sub- 
sidiary conditions, but they have applications in the anal- 
ysis of the P-representation. See Appendix B for the def- 
inition of the P-representation of Gaussian states. 



III. SEPARABILITY AND 
P-REPRESENTATION OF GAUSSIAN STATES 

A. General analysis 

One can bring any given V in (19) by Spj2,R) £g> 
Sp(2,R) transformations to the standard form [4|, [f| (see 



also Appendix A) 



V 



a ci 

a c 2 

ci b 

c 2 b 



(30) 



One may understand (22) and (23) ( and similarly (28) 
and (29)), which hold for any choice of parameters d to 
h, as imposing constraints on the allowed ranges of the 
elements a, 6, c\,ci of the standard form of Vq in (30), 
for example. This interpretation was also adopted by 
Simon [4j. In this interpretation, the full relation (22) 
is more restrictive than the relations (28) and (29) with 
the subsidiary conditions on the parameters d to h. In 
other words, the elements a, 6, Ci,C2 which satisfy (22) 
automatically satisfy (28) and (29), but not the other 
way around. In our analysis below, we interpret these 
relations as constraints on c\, ci for fixed a, b. [23] 

We now analyze the Gaussian states. The existence of 
the well-defined P-representation of two-party Gaussian 
states requires the condition (see Appendix B) 



P = V / > 

2 ~ 



(31) 



namely 



d T Ad + f T Bf + 2d T Cf>-(d T d + f T f) (32) 

for any d and /. By using a special property of the P- 
representation, namely, a special property of the coherent 
state, one can identify V in (21) with the matrix P^ 1 if 
the P-representation exists. See (B13) in Appendix B. 
We thus have 



p- 1 = V = V - -I 
2 



(33) 



or equivalently 

d T Ad + f T Bf + 2d T Cf 

= d T Ad + f T Bf + 2d T Cf - ^(d T d + f T f) (34) 

which is in fact non-negative for any d and /, as is re- 
quired by (32), since V is positive semi-definite by its 
construction. 

We first show that the P-representation condition im- 
plies the separability conditions. When one adds (34) 
to the expression obtained from (34) by the replacement 
of d and / by g and h, respectively, one reproduces the 
separability condition (22) 

d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 
= d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 

+ \{d T d + f T f)+ 1 -(g T g + h T h) 

> d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 



-\d T Jg\ + \f T Jh\, 



(35) 
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where we used 
1 



(d 1 d + g 1 g) + -(f 1 f + h 1 h) 



> 



y/(<Fd)(gTjTjg) + \J (f T f)(h T J T Jh) 

> \d T Jg\ + \f T Jh\. (36) 

This is natural since the P-representation of Gaussian 
states is in fact separable. 

But the converse is not obvious and needs to be 
proved. The separability condition (22) is invariant under 
Sp(2, R)®Sp(2, R) in (25), whereas the P-representation 
condition (31) or (32) is not invariant under Si (g> S 2 £ 
Sp(2, R) <g> Sp(2, R) by noting that 

d T S^S 1 d + f T S 2 r S 2 f + d T d + f T f (37) 

in general. In this sense these two conditions cannot 
be equivalent to each other. In the process of the iden- 
tification (33), the Sp(2,R) ® Sp(2,R) invariance of V 
is lost due to the normal ordering operation in the P- 
representation. See Appendix B. 

The separability condition of Simon Q , which is based 
on the non-negativity of the partially transposed density 
matrix, is given by 



V 



J 
±J 



> 



(38) 



or equivalently (by taking the average of (38) in the form 
v' Mv by the four-component complex vectors v = (d ± 
ig,f±ih)) 

d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 
> \d T Jg\ + \.f T Jh\ (39) 

which is Sp(2,R) ® Sp(2,R) invariant. Since V is non- 
negative, the condition (22) is generally more stringent 
than (39), and thus the condition (38) or (39) is nec- 
essary but not sufficient for separability in general even 
for the case when (22) is sufficient. However, for the P- 
representation of Gaussian states, we have V = P^ 1 = 
V — (1/2)/ as in (33) and thus V is determined by V. 
Moreover at the boundary of the P-representation con- 
dition, the two eigenvalues of V — (1/2)1 after a suitable 
squeezing transformation vanish [22| • For this reason, the 
separability condition (39) can be equivalent to (22) and 
thus can also be sufficient for the separability of Gaus- 
sian states. This is indeed the case as we have shown 
elsewhere by an explicit algebraic analysis [22| . 

There is yet another set of separability conditions for- 
mulated by Duan et al. [f|. To be precise, the condition 
used in Q corresponds to the weaker separability condi- 
tion (29) 



A + JAJ T C ± JCJ T 
C T ±JC T J T B + JBJ T 



> I 



(40) 



applied to the separability condition of Simon (39) by 
setting V = 0. This fact, which has been briefly men- 
tioned elsewhere [22| . is explained in more detail in the 
following. 



B. Quantitative analysis 

We start with the analysis of the P-representation con- 
dition V - [1/2)1 > in (31). We apply this condition 
to the covariancc matrix 



V = S^VoiS- 1 ) 11 

ari ciy/rir 2 

a/ri c 2 /V r i r 2 

ciV r i r 2 br 2 

c 2 / ^jT\r 2 b/r 2 



41) 



which is obtained from the standard form Vq in (30) by 
a transformation S(ri, r 2 ) £ Sp(2, R) ® Sp(2, R) defined 
by 



S(r 1 ,r 2 )S T (r 1 ,r 2 ) = 



1/ri 

n 

l/r 2 

r 2 



(42) 



with suitably chosen n > 1 and r 2 > 1. By evaluating 
eigenvalues of ( 41), one then obtains the conditions for 
V - (1/2)1 > [22J 



(a-^)(6-^)>c?, 
2ri 2r 2 



1 



1 



(a - -n)(b - -r 2 ) > c. 



(43) 



together with 

(«-^) + (*-^)>0, (a 



-r 1 ) + (b--r 2 )>m 



which require a > 1/2 and b > 1/2. The boundary of the 
P-representation conditions (43), namely, the extremal of 
these conditions with respect to r\ and r 2 when regarded 
as the bound on c\ is specified by 



(a- 



2r x 



)(b- 



2r- 



■) 



t: 2 



l(a-7;n)(b--r 2 )} (45) 



together with the constraint which defines the extremal- 
ity [12] 



(an -1/2) _ (br 2 - 1/2) 



(a/n-1/2) {b/n- 1/2)' 
Here we defined 

< t = Ical/lcil < 1 



(46) 



(47) 



by choosing \c 2 \ < \ci\ without loss of generality. The 
conditions (45) and (46) are explicitly solved as [22| 



1 



n 

r 2 



at + b 
1 



a 



bt 



{ab(l-t 2 ) + ^D(a,b,t)}, 

{ab(l - t 2 ) + y/D(a,b,t)} (48) 
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with the auxiliary quantity 



D(a, b, t) = y/a 2 b 2 {l - t 2 ) 2 + t(a + bt){at + b). (49) 

These squeezing parameters are shown to satisfy [22] 

2a > r x > 1, 2b > r 2 > 1 (50) 

for a > 5 and b > \ , to be consistent with (44) . 

For these values of squeezing parameters in (48), we 
have the upper bound on |ci| from (43) by considering 
the square root of (45) [12] 



\ci\ < 



V(ari-l/2)(br 2 -l/2) 

_ V(o/ri-l/2)(6/f2-l/2) 



= ^{[2afe(l + < 2 ) + 1] - 2/DKM)} 1/2 (51) 

The condition (51) defines the P-representation condi- 
tion which characterizes the separable Gaussian states, 
namely, we have the P-representation for any |ci| which 
satisfies (51) for any given a > i, b > h and 1 > t > 0. 

The separability condition (38) or (39) derived by Si- 
mon is written as the algebraic conditions 



A(ab - c\)(ab - c\) > (a 2 + b 2 ) + 2\ Cl c 2 \ 
v / (2a-l)(26-l) > |ci| + |c 2 | 



(52) 



together with a > 1/2 and b > 1/2. The second con- 
dition in (52) is derived from the weaker condition (40) 
applied to the standard form Vq in (30) [H)]; the second 
condition is used to exclude the solutions with |ci| — > oo 
and |ci| — > oo for fixed a and b, which are allowed by the 
first condition in (52) of Simon. These conditions (52) 
are explicitly solved as [22j 



c < 



At 2 



{[2ab{\ + t 2 ) + t\- 2y/D(a,b,t)}, 



c\ < ~{[2ab{l + t 2 ) + 1] - 2y/D{a, b, t)} (53) 

which agree with the P-representation condition (51). 
Any separable state, in particular, the P-representablc 
state satisfies the condition (38) by its construction. The 
relations (51) and (53) now show that any Gaussian 
state which satisfies the condition (38) satisfies the P- 
representation, and thus (38) gives a necessary and suffi- 
cient criterion for the separability of two-party Gaussian 
systems 

We now come to the weaker form of the separability 
condition (40). It is confirmed that the condition (40) 
when applied to the covariance matrix V in (41) gives 
rise to (see also the second condition in (52)) 



Vi + - - l][6ra + - - 1] > v^kil + ^L(54) 



If one imposes the condition (46), the left-hand side of 
(54) becomes 



a b 

[an + 1 ][br 2 + 1] 

n r 2 



(55) 



= v/(an - l/2)(6r 2 - 1/2) + y/(a/n - l/2)(6/r 2 - 1/2). 

See Appendix C. The weaker form of separability condi- 
tion (40) together with (46) thus gives 



vVi - 1/2) (br 2 - 1/2) + y/(a/ ri - l/2)(6/r 2 - 1/2) 
IPal 



> xM^lcil + 



(56) 



which is confirmed to correspond to the separability con- 
dition of Duan et al., eq. (16) in [5j, when converted into 
their notation. 

When one sets |c 2 | = t\ci\ in (56), one obtains the 
condition 



\/(ari 



- l/2)(6ra - 1/2) + v/(a/ri-l/2)(6/r 2 -l/2) 
>[v^+-^=]| Cl |. (57) 



One can confirm that this relation (57) when regarded as 
a bound on |ci| agrees with the P-representation con- 
dition (51) if one uses the first equality in (51). We 
thus conclude that the weaker separability condition (40) 
when combined with the squeezing defined by (48) pro- 
vides a necessary and sufficient criterion for the separa- 
bility of two-party Gaussian states. This proof, which 
was sketched in |22j ]. is much simpler and more explicit 
than the original proof in [f| which was also completed 

in m. 

In the above analysis, we implicitly used the solution 
of (46) 



r 2 (ri) 



46 



[yj (1 — X) 2 + 16& 2 A + (1 - X)] 



(58) 



with X{ri) = (2a /ri — l)/(2ari — 1) which assumes 
X(l) = 1 and X(2a) = 0, and thus r 2 (l) = 1 and 
r 2 (2a) = 2b. 



IV. CONCLUSION 

We have derived a separability criterion for contin- 
uous variable two-party systems without referring to 
the non-negativity of the partially transposed density 
matrix. This condition on the second moment (or 
covariance matrix) is generally more stringent than the 
condition of Simon [4| which is based the non-negativity 
of the partially transposed density matrix. Our criterion 
may thus be useful as a necessary condition for the 
separability of general two-party systems. We have also 
shown that the separability condition used by Duan 
et al. @ is generally weaker than the separability 
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condition of Simon. Nevertheless, at the boundary of 
the P-representation condition for Gaussian two-party 
systems, all these separability conditions become equiv- 
alent to each other if one uses a suitable squeezing 
operation. This fact is most easily shown by using the 
explicit analytic formulas of squeezing parameters which 
establish the equivalence of the Sp(2, R) ® Sp(2, R) in- 
variant separability condition with the P-representation 
condition of Gaussian states [13] ■ 

I thank K. Shiokawa for an informative discussion at 
the very beginning of the present study. 



APPENDIX A: STANDARD FORM OF V 



We recall the elements of Sp(2, R) 



S 



cos V sin t 
— sin 8 cos t 



S 



x 
I 



(Al) 



which satisfy SJS T — J. One can bring V in (17) to the 
standard form 



V 



( a ci 

a c 2 

ci b 

V c 2 b 



(A2) 



by suitable Sp(2,R) ® Sp(2,R) transformations 0, 
real symmetric A and B can be made diagonal by two- 
dimensional rotations with suitable parameters in (Al) 
and then applying the second elements in (Al) with 
suitable parameters x, A and B are made proportional 
to the unit matrix. After these transformations C re- 
mains real. By applying a suitable two-dimensional or- 
thogonal transformation S\ ® S%, which is an element of 
Sp(2, R) ® Sp(2, R), we can diagonalize C 



ci 
c 2 



(A3) 



By this way we arrive at (A2). 



APPENDIX B: P-REPRESENTATION 

We define the generating function of all the correlations 
(or moments) of dynamical variables by 

X(A,7j) = Tr(/3cxp{i(Ai<?i + \ 2 p\ + r\\q 2 + %P2)})(B1) 

where Ai ~ 772 are real numbers. By expanding x(A, 77) 
in powers of Ai ^772, one can generate all the moments 
of dynamical variables. One may generally choose (<7i) = 
(pi) = ((72) = (p2j = 0. Following the convention in this 
field, we define the Gaussian states by 



where V is the covariance matrix in (17), namely, all 
the correlation functions are determined by the second 
moments. Note that c-numbers Ai ~ 772 are commuting 
and thus all the operators are automatically symmetrized 
in (Bl). One can write (Bl) as 

x(A,?7) ^Tr^exp^A^a + Aa 1 + rj*b + r}P)}) (B3) 
with 

a= 77|(<Zi +ipi), b= -^=(q 2 +ip 2 ), 

A= ^=(Ai +7A 2 ), 77= -^=(m (B4) 

The Gaussian state is called P-representable if the den- 
sity matrix is written as 



p = d 2 a d 2 /3P(a,f3)\a,f3)(a,(3\ (B5) 



where \a, (3) is the (over complete) coherent state defined 
by 

a\a, f3) = a\a, 0), b\a, (3) = f3\a, f3) , (a, j3\a, (3) = 1 

(B6) 

or to be explicit 

|a,/3) =e Qat -^ |Q|2 |O)0e /3St -^ l/3|2 |O). (B7) 

Thus the P-representable states are separable. 

By using the density matrix (B5) in (B3) and after 
normal ordering the exponential factor in (B3), we have 

X (A,77) = J d 2 aj d 2 f3P(a,P) 

x exp{i(A*a + Xa* + if (3 + rj/3*)} 
xexp{-i(|A| 2 + |77| 2 )} (B8) 

or, if one combines this expression with (B2) we have 
cx P{-^( A i7 A 2 ,»7i>%)(V r - ^)(Ai, A 2 ,77 1 ,?7 2 ) T } 



d 2 a I d 2 /3P(a,(3) 
x exp{i(A 1 Q: 1 + A 2 a 2 + Vifii + f ?2/32)} 

d 2 !3P{a,(3) 
x exp{?:(A*a + Xa* + r]*(3 + 77/?*)} 



(B9) 



with a = (ai + ia 2 )/V2 and /3 = (/3i + ifo)/y/2. Thus 
P(a, (3) in (B9) is given by 



P(<*,0) 



y/detP 
4tt 2 



(BIO) 



x(A,?7) = exp{--(Ai, A 2 ,77i, 77 2 )V r (Ai, A 2 ,77i, 77 2 ) T }(B2) 



x exp{--(ai, a2,/3 1 ,/3 2 )P(ai, a 2 , (3 X , (3 2 ) T } 
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with a matrix P 



where the equality holds only when the condition 



v-- 2 i>o, 



(Bll) 



which defines the condition for the existence of the well- 
defined P-rcpresentation. 

The P-representation is characterized by the weight 
factor P{a, (3) which is in turn determined by the matrix 
P. We now notice that the right-hand side of the formula 
(B9), when expanded in terms of A's and 77's, generates 
the correlations of the form 

J d 2 a J d 2 f3P(a,l3)a*a (B12) 

= J d 2 a J d 2 f3P(a,(3)(a,P\a' i \a,P)(a,(3\a\a,P), 

for example. This may be compared to (20). By re- 
calling (B5), this (B12) shows that all the second mo- 
ments on the right-hand side of (B9), which are deter- 
mined by P, coincide with V in (20) (if one chooses 
(a 1 *) = (a) = (¥) = (b) = 0). This property establishes 
the special relation (33) of the P-representation, namely, 
V is determined by 



V = P . 

APPENDIX C: PROOF OF (55) 

One can show 

VlOi + n 2 ) - l][(mi + m 2 ) - 1] 



(B13) 



>J[ni- \]\mi -\] + J[n2- \]\m 2 - \] (CI) 



("2 - |) = (m 2 - |) 
(m - i) (mi - I) 



(C2) 



is satisfied. This relation (CI) is established by consid- 
ering 



J'(.v) = \/[m + x(n 2 - ni) - ^][mi + x(m 2 - mi) - i] 



(C3) 



and 



\[{ni - \){m 2 -\)~ («2 - J)K - i)] 2 



x [mi + x(m 2 - mi) - ^] 3/2 
x [ni +x(n 2 -m) - ^] _3/2 < 



(C4) 



except for (C2) for which / (a;) = 0. By using the prop- 
erty of the convex function 2/(1/2) > / (1) + /(0) one es- 
tablishes (CI). If one sets m = ari, n 2 = a/r\, mi = br 2 
and m 2 = b/r 2 , (CI) gives (55) while (C2) corresponds 
to (46). 



[1] A. Einstein, B. Podolsky and N. Rosen, Phys. Rev. 47 
(1935) 777. 

[2] A. Peres, Phys. Rev. Lett. 77 (1996) 1413. 

[3] P. Horodecki, Phys. Lett. A232 (1997) 333. 

[4] R. Simon, Phys. Rev. Lett. 84 (2000) 2726. 

[5] L.M. Duan, G. Giedke, J.I. Cirac and P. Zoller, Phys. 

Rev. Lett. 84 (2000) 2722. 
[6] R.F. Werner and M.M. Wolf, Phys. Rev. Lett.86 (2001) 

3658. 

[7] E. Shchukin and W. Vogel, Phys. Rev. Lett. 95 (2005) 
230502. 

[8] A. Miranowicz and M. Piani, Phys. Rev. Lett. 97 (2006) 
058901. 

[9] B.G. Englert and K. Wodkiewicz, Phys. Rev. A65 (2002) 
054303. 

[10] G. Giedke, B. Kraus, M. Lewenstein, J.I. Cirac, Phys. 

Rev. Lett. 87 (2001) 167904. 
[11] S. Mancini, V. Giovannetti, D. Vitali, P. Tombesi, Phys. 

Rev. Lett. 88 (2002) 120401. 
[12] J. Eisert, S. Scheel, and M.B. Plenio, Phys. Rev. Lett. 89 

(2002) 137903. 



[13] G. Vidal and R.F. Werner, Phys. Rev. A65 (2002) 
032314. 

[14] M.M. Wolf, J. Eisert, and M.B. Plenio, Phys. Rev. 
Lett. 90 (2003) 047904. 

[15] M.G. Raymer, C. Funk, B.C. Sanders, H. de Guise, Phys. 
Rev. A67 (2003) 052104. 

[16] V. Giovannetti, S. Mancini, D. Vitali, P. Tombesi, Phy. 
Rev. A67 (2003) 022320. 

[17] G. Giedke et al., Phys. Rev. Lett. 91 (2003) 107901. 

[18] S. Mancini and S. Severini, Electronic Notes in Theo- 
retical Computer Science 169 (2007) 121, and references 
therein. 

[19] S.L. Braunstein and P. van Loock, Rev. of Mod. Phys. 77 
(2005) 513. 

[20] E.H. Kennard, Z. Phys. 44 (1927) 326. 

[21] HP. Robertson, Phys. Rev. 34 (1929) 163. 

[22] K. Fujikawa, Phys. Rev. A79 (2009) 032334. 

[23] Alternatively, one may understand the relations (22) and 
(23) as an infinite set of uncertainty relations (and their 
variants) for each given set of parameters d to h. 



